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NATIONAL. ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 


A STUDY OF THE MOTION AND AERODYNAMIC HEATING 
OF MISSILES ENTERING THE EARTH'S ATMOSPHERE 
AT HIGH SUPERSONIC SPEEDS 


By H. Julian Allen and A. J. Eggers, Jr. 
SUMMARY 


A simplified analysis is made of the velocity and deceleration 
history of missiles entering the earth's atmosphere at high supersonic 
speeds. It is found that, in general, the gravity force is negligible 
compared to the aerodynamic. drag force and, hence, that the trajectory 
is essentially a straight line. A constant drag coefficient and an 
exponential variation of density with altitude are assumed and general- 
ized curves for the variation of missile speed and deceleration with 
altitude are obtained. A curious finding is that the maximum decelera- 
tion is independent of physical characteristics of a missile (e.g., mass, 
size, and drag coefficient) and is determined only by entry speed and 
Plight-path angle, provided this deceleration occurs before impact. 

This provision is satisfied by missiles presently of more usual interest. 


The results of the motion analysis are employed to determine means 
available to the designer for minimizing serodynamic heating. Emphasis 
is placed upon the convective-heating problem including not only the 
total heat transfer but also the maximum average and local rates of 
heat transfer per unit area. It is found that if a missile is so heavy 
as to be retarded only slightly by aerodynamic drag, irrespective of 
the magnitude of the drag force, then convective heating is minimized 
by minimizing the total shear force acting on the body. This condition 
is achieved by employing shapes with a low pressure drag. On the other 
hand, if a missile is so Light as to be decelerated to relatively low 
speeds, even if acted upon by low drag forces, then convective heating 
is minimized by employing shapes with a high pressure drag, thereby 
maximizing the amount of heat delivered to the atmosphere and minimizing 
the amount delivered to the body in the deceleration process. Blunt 
shapes appear superior to slender shapes from the standpoint of having 
lower maximum convective heat-transfer rates in the region of the nose. 
The maximum average heat-transfer rate per unit area can be reduced by 
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employing either slender or blunt shapes rather than shapes of inter- 
mediate slenderness. Generally, the blunt shape with high pressure 
drag would appear to offer considerable promise of minimizing the heat 
transfer to missiles of the sizes, weights, and speeds presently of 
interest. 


INTRODUCTION 


In the design of long-range rocket missiles of the ballistic type, 
one of the most difficult phases of flight the designer must cope with 
is the re-entry into the earth's atmosphere, wherein the aerodynamic 
heating associated with the high flight speeds of such missiles is 
intense. The air temperature in the boundary layer may reach values in 
the tens of thousands of degrees Fahrenheit which, combined with the 
high surface shear, promotes very great convective heat transfer to the 
surface. Heat-absorbent material must therefore be provided to prevent 
destruction of the essential elements of the missile. It is a charac- 
teristic of long-range rockets that for every pound of material which 
is carried to burn-out,” many pounds of fuel are required in the 
booster to obtain the flight range. It is clear, therefore, that the 
amount of material added to protect the warhead from excessive aero- 
dynamic heating must be minimized in order to keep the take-off weight 
to a practicable value. The importance of reducing the heat transferred 
to the missile to the least amount is thus evident. 


For missiles designed to absorb the heat within the solid surface 
of the missile shell, a factor which may be important, in addition to 
the total amount of heat transferred, is the rate at which it is trans- 
ferred since there is a maximum rate at which the surface material can 
safely conduct the heat within itself. An excessively high time rate 
of heat input may promote such Large temperature differences as to 
cause spalling of the surface, and thus result in loss of valuable heat- 
absorbent material, or even structural failure as a result of stresses 
induced by the temperature gradients. 


For missiles designed to absorb the heat with liquid coolants 
(e.g., by "sweat cooling" where the surface heat-transfer rate is high, 
or by circulating liquid coolants within the shell where the surface 
heat-transfer rate is lower), the time rate of heat transfer is simi- 
larly of interest since it determines the required liquid pumping rate. 


These heating problems, of course, have been given considerable 
study in connection with the design of particular missiles, but these 
studies are very detailed in scope. There has been need for a general- 
tized heating analysis intended to show in the broad sense the means 
available for minimizing the heating problems. Wagner, reference l, 
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made a step toward satisfying this need by developing a laudably simple 
motion analysis. This analysis was not generalized, however, since it 
was his purpose to study the motion and heating of a particular missile. 


Et is the purpose of this report to simplify and generalize the 
analysis of the heating problem in order that the salient features of 
this problem will be made clear so that successful solutions of the 
problem will suggest themselves, 


A motion analysis, having the basic character of Wagner's approach, 
precedes the heating analysis. The generalized results of this analysis 
are of considerable interest in themselves and, accordingly, are treated 
in detail. 


ANALYSIS 
Motion of the Body 
Consider a body of mass m entering the atmosphere from great 
height. If, at any altitude y, the speed is V and the angle of 


approach is 8 to the horizontal (see sketch), the parametric equa- 
tions of motion can be written 


Impact point (O, O) 


*Properly, the analysis should consider those effects resulting from the 
fact that the earth is a rotating sphere, but since the altitude range 
for which drag effects are important is less than 1 percent of the 
radius of the earth, the rectilinear treatment given in this analysis 
is permissible. 
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where 


Cp drag coefficient, dimensionless 


<j 


speed, ft/sec 

A reference area for drag evaluation, sq ft 
m mass of the body, slugs 

p mass density of the air, slugs/ft? 

g acceleration of gravity, ft/sec? 


X,y horizontal and vertical distance from the point of impact with 
the earth, ft 


9 angle between the flight path and the horizontal, deg 

(See Appendix A for complete list of symbols.) 

In general, the drag coefficient varies with Mach number and 
Reynolăs number, while the density and, to a very minor extent, the 
acceleratlon of gravity vary with altitude. Hence it 1s clear that 


exact solution of these equations 1s formidable. Let us first, then, 
consider the following simplified case: 


1. The body descends vertically. 
2, The drag coefficient is constant.” 
3. The acceleration of gravity is constant.“ 


4. The density as a function of altitude is given by the relation 
p = pe PY (2) 


where P, and fp are constants* 


“As is well known, this assumption is generally of good accuracy at the 
high Mach numbers under consideration, at least as long as the total 
drag is largely pressure drag. 

SThe acceleration of gravity decreases by only 1 percent for every 
100,000-foot increase in altitude. 

*This relation is consistent with the assumption of an isothermal 
atmosphere. - 


ie. 
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Equations (1) then reduce to the single equation 


ay dv CDPoA -By 
=- œ R L US 
ate at Bee a e (3) 
Noting that 

AV... yY 

at dy 
we let 

Z = VS 


ana equation (3) becomes the linear differential equation 


C 
dy m 


which has the well-known solution 
CP. A C PA 
Do - D = 
Pi mm” e PY ay -f a e PYay 
A= e - 2g i e 


Performing the integrations, we obtain as the solution of this relation 


dy + const, 


n 
_SpPoA -ayl < Ga "By 
7, = V2=e fe as E ee ØRE Sey + const. (5) 


nia 


so that the deceleration becomes, in terms of gravity acceleration, 


av Cpe A CpPoâ J 
tdi = - By co -By 
Gg SDP of „Bye Bm = 2g 3 NB =- Sgytconst.| -1 
g 2mg B nja | 
s1 (6) 
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As an example, consider the vertical descent of a solid iron sphere 
having a diameter of 1 foot. For a sphere the drag coefficient may be 
taken as unity, based on the frontal area for all Mach numbers greater 
than about 1.4. In equation (2), which describes the variation of 
density with altitude, the constants should clearly be so chosen as to 
give accurate values pf the density over the range of altitudes for 
which the deceleration is large. It is seen in figure 1 that for 


po = 0.0034 slugs/ft® 


and 
B = AL. fl 
22,000 
Which yields 
eee 2 
p = 0.0034 e 22,000 (7) 


the calculated density is in good agreement with the NACA standard 
atmosphere values obtained from references 2 and 3 for the altitude 
range from 20,000 to 180,000 feet. These relations have been used in 
calculating the velocity and deceleration of the sphere for various 
altitudes, assuming vertical entrance velocities of 10,000, 20,000, 

and 30,000 feet per second at 40 miles altitude which, for these cases, 
may be considered the “outer reach" of the atmosphere. The results of 
these calculations are presented as the solid curves in figures 2 and 3. 


It is seen in figure 3 that for the high entrance speeds considered, 
the decelerations reach large values compared to the acceleration of 
gravity. This suggests that the gravity term in equation (3) may be 
neglected without seriously affecting the results.5 When this term is 
neglected the equation of motion becomes 


C 
_ I y WW _ SDA TPY ye (8) 
dt dy 


“It is usual to neglect the gravity acceleration a priorl (see e.g., 


refs. 1 and k.) 
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Integration gives 


or 


V = const. X e 2pm 
At the altitude of 40 miles it can readily be shown that the term 


ȘI CDPoA a BY 
e 2pm 


is very nearly unity so that the velocity may be written 


OPPA -By 
Y < Ve 2pm (9) 


and 


av CP. A 
BE Po? gy — Se PY 
ee a (10) 


where VE is the entrance speed. 


By use of equations (9) and (10) the vertical-decent speeds and 
decelerations for the 1-foot-diameter sphere previously considered have 
been calculated for the same entrance speeds. The results are shown as 
the dashed curves in figures 2 and 3. It is seen that these approximate 
calculations agree very well with those based on the more complete 
equation of motion (eq. (3)). 


The above finding is important, for it indicates that in the gen- 


eral case, wherein the body enters the atmosphere at high speed at 
angle 0p to the horizontal, the gravity term, provided Om is not 


AE 
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too small, may be neglected in equation (1) to yield 


a2y  CphPVTA sin 6, 
att Om 
(11) 


dx _ CppVEA cos On 
at? om 


so that the flight path is essentially a straight line (1.e., 6 = Op), 
and the resultant deceleration equation becomes 


ve 
_ăv oA (12) 


Now, again, if the density relation given by equation (2) is used and 
1t is noted that 


_ 1 
V= dt or -N aysi 6, A 
sin 0p dt 


av CDPo 2” BY 


U z 2m sin On dy 
which can be integrated to yleid 
D 
age (23) 
V = VE € SASAE 
and the deceleration is then 
av 2 _ CpPoA e PY 
LG. Ka SE eNe ME (14) 
E omg 


[map prea 
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The altitude yı at which the maximum deceleration occurs is found 
from this relation to be 


Cp? of 


1 
= = m =m 15 
Yi B am sin Or ( ) 
If yı is positive, the velocity Va (from eqs. (13) and (15)) at 
which the maximum deceleration occurs becomes 
Va = Vg 3 = 0.61 Vg (16) 


and the value of the maximum deceleration is 


GV ay 2 
_{ dt _ NEN, NG sin % (17) 
8 g 2ge 
MAX 1 


If equations (13) and (14) are rewritten to make the altitude 
reference point yi rather than zero, then 


CpPoh _.-Blyi+Ay) 


v= Ve 2pm sin 0p 
and 
av = _ _SpPoA _ -B(Yz+Ay) 
at _ PoE. -B(y1+4y) „ Pm sin 65 
E omg 


respectively, where Ay is the change in altitude from yi. Substitu- 
tion of equation (15) into these expressions can readily be shown to 
give 


Z re” = F' (BAY) (18) 
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and 


av/at 
(Se) = e” PAY „(ae AY) = F"( pay) 
(=) 
& 1 


Equations (18) ana (19) are generalized expressions for velocity and 
deceleration for bodies of constant drag coefficient and, together with 
equations (15) and (17), can be used to determine the variation of 
these quantities with altitude for specific cases. The dependence of 
F'(PAy) and F"(BAy) on My is shown in figure 1. 


(19) 


The maximum deceleration and the velocity for maximum decelera- 
tion as given by equations (17) and (16) apply only if the altitude yu, 
given by equation (15), is positive. Otherwise the maximum decelera- 
tion in flight occurs at sea level with the velocity (see eq. (13)) 


papal 
V = Vo = Ye e S sin Ôg (20) 
and has the value 
ay. av Zi we oe 
Ft) (Æ | . Doe - Rein og (21) 
6 g eng 


Heating of the Body 


It was noted previously that for practicable rocket missiles, it 
is vital that the weight of the missile be kept to a minimum. The 
total heat transferred to a missile from the air must be absorbed by 
some "coolant" material. Since this material has a maximum allowable 
temperature, it follows that 1t can accept only a given amount of heat 
per unit weight. Hence, the total heat input to the missile must be 
kept at a minimun for minimum missile weight. 
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Often the coolant material is simply the shell of the missile and 
as such must provide the structural strength and rigidity for the 
missile as well. The strength of the structure is dictated, in part, 
by the stresses induced by temperature gradients within the shell. 
Since these temperature gradients are proportional to the time rate of 
heat input, the maximum time rate of heat input is important in missile 
design. The heating, of course, varies along the surface but, since 
the shell transmits heat along as well as through itself, the strength 
of the structure as a whole may be determined by the maximum value of 
the average heat-transfer rate over the surface. This is simply the 
maximum value of the time rate of heat input per unit area. On the 
other hand, the structural strength at local points on the surface may 
be determined primarily by the local rate of heat imput. Hence, the 
maximum time rate of heat input per unit area at the surface element 
where the heat transfer is greatest may also be of importance in design. 


If liquid cooling is employed, the maximum surface heat-transfer 
rates retain their significance but, now, in the sense that they dic- 
tate such requirements as maximum coolant pumping rate, or perhaps 
shell porosity as well in the case of sweat cooling. Whichever the 
case, in the analysis to follow, these elements of the heating problem 
will be treated: 


i. The total heat input 
2. The maximum time rate of average heat input per unit area 
3. The maximum time rate of local heat input per unit area 


Since it is the primary function of this report to study means 
available to the missile designer to minimize the heating problem, the 
analysis is simplified to facilitate comparison of the relative heating 
of one missile with respect to another - accurate determination of the 
absolute heating of individual missiles is not attempted. With this 
point in mind, the following assumptions, discussed in Appendix B, are 
made: 


l. Convective heat transfer predominates (1.e., radiation effects 
are negligible). 

2, Effects of gaseous imperfections may be neglected. 

3. Shock-wave boundary-layer interaction may be neglected. 

h. Reynolds! analogy is applicable. 

5. The Prandtl number is unity. 


Total heat input .- The time rate of convective heat transfer from 


the air to any element of surface of the body may be expressed by the 
well-known relation 


SE = by (Ir - Ty); (22) 


ETT ng 
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where 
H heat transferred per unit area, ft-1b/ft* 


h convective heat-transfer coefficient, ee 
ft? sec OR 
T, recovery temperature, OR 

Ty temperature of the wall, OR 

t time, sec 


and the subscript 2 denotes local conditions at any element of the 
surface dB. l 


It-is convenient in part of this analysis to determine the heating 
as a function of altitude. To this end, noting that 


åt = ae 
V sin Op 


we see that equation (22) may be written 


bh-(Tr - Ty) 
dy V sin êp 


With the assumption that the Prandtl number is unity, the recovery tem- 


perature is 
Yel o Vol 
Tp Pa (14 TF m) = (a + 2 w) 


where 
M Mach number at the altitude y, dimensionless 


y the ratio of specific heat at constant pressure to that at constant 
volume, Cy/Cys dimensionless 


T static temperature at the altitude y, OR 
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so that 


-1 
(Te - Ty), = T- Ty, tM 


It is seen that for large values of the Mach number, which is the case 
of principal interest, the third term is large compared to reasonably 
allowable values of T - Ty. It will therefore be assumed that T - Ty 
is negligible so that 


(Tr - Tu), oe sg (24) 
Moreover, since 
2 
Yn lia 
(7-1)C, 
equation (2) may be written 
y2 ) 
T T se 2 
(Fr wh 2C (25 


Now the local heat-transfer coefficient h, is, by Reynolds' analogy, 
for the assumed Prandtl number of unity 


Ra = 3 Ce, Cp, Pz Vy (26) 


Where Cr. is the local skin-friction coefficient based on conditions 


P13 Vz, etc., just outside the boundary layer. Thus, since (To - Ty) 
is essentially constant over the entire surface S, the rate of total 


STt should be noted that without this assumption, the heat-input deter- 
mination would be greatly complicated since the changing wall tem- 
perature with altitude would have to be considered to obtain the heat 
input (see e.g., ref. 1). For high-speed missiles which maintain 
high speed during descent, the assumption 1s obviously permissible. 
Even for high-speed missiles which finally decelerate to low speeds, 
the assumption is generally still adequate since the total heat input 
is largely determined by the heat transfer during the high-speed 
portion of flight. 


CG au 
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heat transfer with altitude becomes from equations (23) through (26) 


dQ f2 Y 
= | 24S = - G P; V; ds 
dy op in Og | fy CP, Pr Vi 


where Q is the heat transferred to the whole surface 5. This equa- 
tion may be written 


dy k sin CM 


wherein Cp is set equal to Cp and 
L 


3 [ e, (Ge m 
S 


The parameter Cr' is termed "the equivalent friction coefficient," 
and will be assumed constant,” independent of altitude, again on the 
premise that relative rather than absolute heating is of interest. 
With equations (2) and (13), then, equation (27) is written 


CpPoA -PY 
2 NATA. = 
E or Ng PA (29) 


Comparison of equation (29) with equation (14) shows that the 
altitude rate of heat transfer 1s directly proportional to the 


rrr 

7This assumption would appear poor at first glance since the Mach number 
and Reynolds number variations are so large. Analysis has indicated, 
however, that the effects of Mach number and Reynolds number variation 
are nearly compensating. The variation in Cr’ for typical conical 
missiles was found to be, at most, about 50 percent from the maximum 
Cp! in the altitude range in which 80 percent of the heat is 
transferred. | 

kila 
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deceleration, so that 


_aQ/ay_ (E ) (30) 


eS ma 2 sin Or sin Op 


and therefore the maximun altitude rate of heat transfer occurs at the 
altitude yı (see eq. (15)) and is given by 


(2 a) - (8) = - E (=) (31) 


It follows, of course, that the altitude rate of heat transfer varies 
with incremental change in altitude from yi in the same manner as 
deceleration, and thus (see eq. (19)) 


(aa/ay)ay 


(aa/ay), F (PAY) (32) 


The total heat input to the body at impact follows from equation (29) 
(integrating over the limits O< y< œ) and is 


__CpPoA _ 
C£'S - Bm sin Bm sin Om 
T a Sa =) nV,2\1 - e œ (33) 
The impact velocity, V, (the velocity of body at y = 0), is 


Dron 


— 2fm sin Oy, 


so that equation (33) may be written in the alternative form 


a= (SE) cae - ve?) (34) 


| mil 
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Maximum time rate of average heat input per unit area.- To deter- 
mine the time rate of average heat transfer per unit area, equations (25), 
(26), and (28) with equation (22) may be shown to give 


diay 
dt 


a 


= m Cp! py? (35) 


= 


Which, together with equations (2) and (13), becomes at altitude y 


Cet Vane <E e 
= 2 
dHay  “FTo'E aY e SIN Dn (36) 


dt ), 


The maximun time rate of average heat transfer per unit area is found 
from this expression to be 


Tar) ar T = E (= =) ma” sin 0g (37) 


max 


and it occurs at the altitude 


DE eaa (38) 
where the velocity is 
a 
Va = Vpe 3 = 0.72 Vp (39) 


As with altitude rate of heat transfer, it can be shown that 


(dHay/dt),, 


= N 40 
(4Ezy/dt),, pama ES 


ma, 
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Equations (37), (38), and (39) apply if the altitude for maximum time 
rate of average heat transfer per unit area occurs above sea level. 
If yo, by equation (38), is negative, then this rate occurs at sea 
level and is, from equation (36), 


oP oh 
day ) GHaxy 3 25 C TVE pa = pan sin On ( 41) 
dt max dt œ h 


Maximum time rate of local heat input per unit area.- The elemen- 
tal surface which is subject to the greatest heat transfer per unit 


are& is, except in unusual cases, the tip of the missile nose which 
first meets the air. It seems unlikely that a pointed nose will be of 
practical interest for high-speed missiles since not only is the local 
heat-transfer rate exceedingly large in this case, but the capacity 
for heat retention is small. Thus a truly pointed nose would burn 
away. Body shapes of interest for high-speed missiles would more prob- 
ably, then, be those with nose shapes having nearly hemispherical tips. 
The following analysis applies at such tips. 


It is well known that for any truly blunt body, the bow shock wave 
is detached and there exists a stagnation point at the nose. Consider 
conditions at this point and assume that the local radius of curvature 
of the body is o (see sketch). 
The bow shock wave is normal to 
the stagnation streamline and 
converts the supersonic flow 
ahead of the shock to a low 
subsonic speed flow at high 
static temperature downstream 
of the shock. Thus, it is sug- 
gested that conditions near the 
stagnation point may be investi- 
gated by treating the nose sec- 
tion as 1f it were a segment of 
a sphere in a subsonic flow field. 


Bow shock wave 


Missile 
nose 
Stagnation 


streamline 


The heat-transfer rate per 
unit area at the stagnation point 
is given by the relation 


di, Nur ky(Ty - Tr) ko 
"aka aaa (42) 


18 ÎI Bhi, NACA RM A53D28 


where kr is the thermal conductivity of the gas at the recovery tem- 
perature (i.e., total temperature) Ty, and Nu, is the Nusselt number 
ot the flow. If the flow is assumed to be laminar and incompressible,8 
Nu, is given, according to reference 5, by the relationship 


2 £g 
Nur = 0.934 Reg? Pr 5 


We retain the assumption that the Prandtl number is unity, note that 
Re, = PVo/h,., and substitute equation (25) into equation (42) to obtain 


Bso] —— 


reg = (43) 


Now 1t is well known that at the high temperatures of interest here, 
the coefficient of viscosity WU, varies nearly as the square root of 
the absolute temperature and is given by the relation 


i 
Hy = 2.31 x 107 Z 


If this expression is combined with equation (25) (neglecting Ty), 
equation (43) may then be written? 


Ze = 6.8 x 1078 = US (h) 


Sme assumption of constant density certainly may invalidate this anal- 
ysis for any quantitative study of the relatively "cold-vell" flows 
of interest here. For the purpose of studying relative heat transfer 
it should, however, prove adequate. 

Had wall conditions rather than recovery conditions been used in the 
development of equation (44), the relation 


dis -7 JP 5 
rai = 1.1 x 10 Jex 


would have been obtained assuming a linear variation of viscosity with 
temperature (to be consistent with the assumption of a cool wall). 
This relation would give somewhat higher heat-transfer rates per unit 
area than equation (44) at velocities greater than about 3600 feet per 
second. 


ED © 


Be ae 
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which, when combined with equations (2) ana (13), becomes 


m SCnD-A - By 
dH -6 — sw 2 sin 8 
ZE = 6.8x 10) Lue 2 e i fg (45) 
The maximum value of dH,/dt can readily be shown to be 
dH dH sin 
= = = = 6.8 x 1076 po a Ce a (46) 
max 3 3eCpcA 
Which occurs at the altitude 
3Cr-PgA 
yy = bm ( RA) (47) 
corresponding to the velocity 
— + 
Va = Vp e 82 0.85 Ve (1:8) 


The manner in which the heat-transfer rate per unit area at the stagna- 
tion point varies with incremental change in altitude from y, can be 
shown to be | 


(45/46), Me md Z [nd 
e e 2 e = prt BAY lL 
(dH, /ăât), RI (49) 


The dependence of F"!'(RAy) on My is shown in figure |. 

Equation (16) applies only if Yg is above sea level. If ya 
from equation (47), is negative, then the maximum heat-trensfer rate 
per unit area at the stagnation point occurs at sea level and is 


SCDPo4 
2) = -8 / Po 3 7 2pm sin 8 


ON TA 
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DISCUSSION 


Motion 


The motion study shows some important features about the high-speed 
Scent of missiles through the atmosphere. The major assumptions of 
this analysis were that the drag coefficient was constant and the den- 
sity varied exponentially with altitude. It was found that the decel- 
eration due to drag was generally large compared to the acceleration 
of gravity and, consequently, that the acceleration of gravity could be 
neglected in the differential equations of motion. The flight path was 
then seen to be a straight line, the missile maintaining the flight- 
path angle it had at entry to the atmosphere. 


For missiles presently of more usual interest, the maximum decel- 
eration occurs at altitude. One of the most interesting features of 
the flight of such a missile is that the maximum deceleration is inde- 
pendent of physical characteristics (such as mass, size, and drag coef- 
ficient of the missile), being dependent only on the entry speed and 
flight-path angle (see eq. (17)). The missile speed at maximum decel- 
eration (eq. (16)) bears a fixed relation to the entrance speed (61 per- 
cent of entrance speed), while the corresponding altitude (eq. (15)) 
depends on the physical characteristics and the flight-path angle but 
not on the entrance speed. It is also notable that for a given incre- 
mental change in altitude from the altitude for maximm deceleration, 
the deceleration and speed bear fixed ratios to the maximum decelera- 
tion and the entry speed, respectively (see Pig. 4 and eqs. (19) 
and (18)), hence, the deceleration and speed variation with altitude 
can readily be determined. 


If the missile is very heavy, the calculated altitude for maximum 
deceleration (eq. (15)) may be fictitious (i.e., this altitude is neg- 
ative) so the maximum deceleration in flight, which occurs just before 
impact at sea level, is less than that calculated by equation (17) and 
is dependent on the body characteristics as well as the entry speed 
and flight-path angle (see eq. (21)). However, the variation of speed 
and deceleration with altitude from the fictitious altitude given by 
equation (15) can still be obtained from figure 1. 


Heating 


Total heat input.- In the heating analysis, a number of simplify- 
ing assumptions were made which should limit its applicability to the 
determination of relative values of heating at hypersonic speeds. It 
is in this relative sense that the following discussion pertains. 
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In considering the total heat transferred by convection to a 
missile, it is evident from equation (33) that the course the designer 
should take to obtain the least heating is affected by the value of 
the factor 


CpP,A 


= B (51) 
pm sin On 


To illustrate, first consider the case of a "relatively heavy" 
missile for which this factor is small compared to unity (the term 
"relatively heavy” is used to denote that the denominator involving 
the mass is very large as compared to the numerator involving the drag 
per unit dynamic pressure, CpA). Then 


map CA 
a pm sin 6p 


is small compared to 1. If this function is expanded in series and 
only the leading term retained, equation (33) becomes 


'ap 2 
„~ CE ' BP, VE 


= (52) 
4B sin Op 


For the relatively heavy missile, then, the least heat will be trans- 
ferred when Cp'S is a minimum -that is to say, when the total shear 
force acting on the body is a minimum. This result is as would be 
expected, if one notes that requiring B<<1 is tantamount to requir- 
ing the missile to be so heavy that it is retarded only slightly by 
aerodynamic drag in its motion through the atmosphere. Hence, the heat 
input to the missile is simply proportional to the shear force. 


Now let us consider the case when B>>1, or, in other words, 
when this missile is “relatively light." In this event, 


CpPoA 
= sin @ 
l-e în E Z1 


and equation (31) can be approximated 


CP'S 


TEA (E) (53) 
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For the relatively light missile, then, the least convective heating is 
obtained when Cr'S/CpA is a minimum. This is at first glance a 
rather surprising result, for it indicates that the heating is reduced 
by increasing the total drag, provided the equivalent frictional drag 
is not increased proportionately as fast. Physically, this anomoly is 
resolved if the problem is viewed in the following way: The missile 
entering the atmosphere has the kinetic energy 1/2 mv,” but, if 


. fm sin Om B 
VE 


is small, then nearly all its entrance kinetic energy is lost, due to 
the action of aerodynamic forces, and must appear as heating of both 
the atmosphere and the missile. The fraction of the total heat which 
is given to the missile is,?° from equation (33), 


a (=) 


Thus, by keeping this ratio a minimum, as much as possible of the 
energy is given to the atmosphere and the missile heating is therefore 
least. i DU 


In order to illustrate these considerations in greater detail, 
calculations have been made using the previously developed equations 
to determine the heat transferred by convection to a series of conical 
missiles. Two classes of missiles have been considered. Missiles in 
the first class were required to have a base area of 10 square feet. 
Missiles in the second class were required to have a volume of 16 cubic 
feet. Gross weights of O, 1,000, 5,000, 10,000, and œ pounds have 
been assumed, and the entrance angle, @,, has been taken as 30° of arc 
in all cases. Missile heating, up to the time of impact, has been 
calculated as a function of cone angle for entrance speeds of 10,000, 
20,000, and 30,000 feet per second. In these calculations the pressure 
drag coefficient was taken as constant for a particular cone at the 
value corresponding to the entrance Mach number (a value of Tp = 500° R 
was assumed throughout). These coefficients were determined from ref- 
erence 6 for cone angles of 10° and greater. For cone angles less 
than 10°, reference 7 was employed to determine these coefficients 


tONote that even if all the drag is frictional drag, only half the heat 
is transferred to the body. The other half is contained in the 
boundary layer and is left in the air in the body wake. 


NACA RM A53D28 AED 23 


(base drag was neglected in all cases). The total drag coefficient 
was taken as the sum o? the pressure drag coefficient plus the skin- 
friction coefficient, the latter coefficient being taken at its value 
for maximum total heat-input rate with altitude. The boundary layer 
was assumed to be wholly turbulent since the Reynolds number, based on 
length of run along the surface of a cone and local conditions just 
outside the boundary layer, was always greater than about 6 x 10° and, 
in fact,was of the order of billions for the more slender cones. 
Turbulent-boundary~layer data were obtained from references 8 and 9, 
and Sutherland's law for the variation of viscosity with temperature 
was used in obtaining "equivalent flat-plate" heat-transfer coefficients. 


Missile heating calculated in this manner for the fixed-base-area 
and fixed-volume cones is presented in figures 5 and 6, respectively. 
Curves for missiles having densities greater than steel are considered 
improbable and are shown as dashed lines. It is clear that for both 
classes of bodies, when the missile is relatively heavy, the optimm 
solution is obtained by making Cr'S as small as possible (small cone 
angle case) and this optimm is accentuated with increase in speed. On 
the other hand, when the missile is relatively light, reduced heating 
is obtained by making Ce'S/CpA as small as possible (the large cone 
angle case). It is noted also that, in general, the advantage of 
reduced heating of the relatively light, blunt cones is more pronounced 
in the fixed-base-area case than in the fixed-volume case. 


Maximum time rate of average heat input per unit area.- It was 
previously noted that.the maximum time rate of average heat input per 


unit area may be of serious importance in determining the structural 
integrity of missiles entering the atmosphere at high speeds.'31 In 
order to illustrate this fact, consider the case of a missile having a 
shell made of solid material and assume that the rate of heat transfer 
per unit area does not vary rapidly from one surface element to the 
next. Then the rate of transfer of heat along the shell will be small 
compared with the rate of transfer through the shell. The shell stress 
due to heat transfer is that resulting from the tendency toward differ- 
ential expansion through the shell and it is proportional to dI fan 
where T is the temperature at any point n within the shell and ņ 
is measured perpendicular from the shell surface. We define En as 
the thermal conductivity of the shell material; then the rate at which 
heat transfers through the shell per unit area is (az, /an) and this 
must, at 1 = O, equal the rate of heat input per unit surface area. 
For the missile considered as a whole, the maximum value of the average 
thermal stress in the shell is a measure of the over-all structural 


11mmis is the common case when the shell material acts as structural 
support and must also transport or absorb the heat. 


= an 
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integrity and the maximum value of this stress will occur at the sur- 
face when 


is a maximum. 


The course the designer should take to minimize the thermal stress 
for the missile as a whole is dependent, as for the case of total heat 
input, upon whether the missile is relatively heavy or light. For the 
relatively heavy missile the value of B, given by equation (51), is 
small compared to unity. The maximum value of the average thermal 
stress in this case is proportional to (see eq. (41)) 


Per ) mo Ce! Po Vp? 
ara” akan aaa (54) 


and, hence, the least average thermal stress is obtained by making Cp! 
a minimum. On the other hand, for the relatively light missile the 
maximum value of the average thermal stress is proportional to (see 


eq. (37)) 
Tar ) = ( 2 Pu ea VE 2 (55) 
do, CHA 6e | 


and, hence, the least average thermal stress occurs when Cr’ [Cpa 
is a minimum. 


In order to illustrate these considerations in greater detail, 
the maximum values of the time rate of average heat input per unit area 
have been calculated for the constant-base-area and the constant-volume 
cones previously discussed in the section on total heat input. These 
values were determined in much the same manner as those of total heat 
input, with the exception that Cp! was evaluated at y, (rather 
than ND? given by equation (38) when it applies, and otherwise at 

. The results are shown in figures 7 and 8. It is seen that 

Be Aiurea values of average thermal stress are reduced for both the 
slender cones and blunt cones as compared to the relatively large 
values of this stress experienced by cones of intermediate slenderness. 


Maximum time rate of local heat input per unit area.- Perhaps even 


more important than the maximum value of the average shell stress is 
the maximum stress that occurs in the shell at the surface element of 
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the missile nose,? where the local heat-transfer rate is probably the 
greatest, for, in general, this latter stress is many times larger. 

In fact, this rate of local heat input can be so large as to promote 
temperature gradients through the shell that are intolerable even with 
the most highly conductive materials (copper, silver, etc.).7° Thus 
some additional means of cooling, such as sweat cooling, may, in any 
case, be required in this region. 


It was stated previously that pointed-nose bodies are undesirable 
due, in part, to the fact that the local heat-transfer rate per unit 
area. at the tip is excessive. The validity of this statement is 
demonstrated by the results of the analysis. It is clear (see eq. (!Hr)) 
that since the local transfer rate varies inversely with the square 
root of the tip radius, not only should pointed bodies be avoided, but 
the rounded nose should have as large a radius as possible. The ques- 
tion then arises; if the nose radius is arbitrarily Tixed, what course 
is available to the missile designer to minimize the problem of Local 
heating at the stagnation point? From both equations (46) and (50), 
it is seen that for an arbitrary nose radius, if the mass, entry speed, 
and flight-path angle are fixed, then the only way to reduce the stag- 
nation rate of heat input per unit area is to increase the product Cya. 
In fact, a relative stagnation-point heat-transfer rate per unit area, Y, 
may be expressed in terms of B (see eq. (51)), if it is defined as the 
ratio of the maximum stagnation-point heat-transfer rate per unit area 
for a given missile to the maximum rate the same missile would experi- 
ence if it were infinitely heavy. For the infinitely heavy missile, 
the maximum rate occurs at sea level and is (see eq. (50)) 


p 
-6 ges) 3 


so that from equation (50) 


-amemo à -ÊB 
¥ =e = e (56) 


“In this report we are concerned only with bodies. If wings or sta- 
bilizers are used, their leading edges are similarly surface elements 
which experience intense heat transfer. The heating problem with 
wings and stabilizers is, in fact, so serious at very high speeds 
that their use as lifting surfaces appears, at present, inadvisable. 

1Ssee reference 1 for further discussion. 


Gites. 
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Tu 


if the given missile also attains its maximun rate at sea level 
(i.e., Ya = 03 eq. (47)): whereas 


pm sin Om 1 i 
¥ = aaa aaa aa - 2 
/ POP LA => (57) 


if the given missile attains its maximum rate above sea level (eq. (16), 
Ya positive). The variation of y with 1/B is shown in figure 9. 
Clearly, the high pressure drag shape has the advantage over the slender 
shape in this respect. 


In order to illustrate these considerations in greater detail, 
again consider the constant-base-area and congtant-volume cones dis- 
cussed earlier. Assume the pointed tips of all the cones are replaced 
by spherical tips of the same radius o. The relative effect of vary- 
ing the cone angle on the stagnation-point heating can then be assessed 
by determining the variation of the product 


NG 7) 
dt /max 


This product has been calculated for the various cones, assuming Cp 

to be unaffected by the addition of the hemispherical tip (the tip i 
radius may be arbitrarily small), and the results are shown in figures 10 

and 11. It is seen again that the missiles having large cone angle 

(high drag coefficient) are considerably superior. = 


DESIGN CONSIDERATIONS AND CONCLUDING REMARKS 


In the foregoing analysis and discussion, two aspects of the heat- 
ing problem for missiles entering the atmosphere were treated. The i | B 
first concerned the total heat absorbed by the missile and was related 
to the coolant required to prevent its disintegration. It was found 
that if a missile were relatively light, the least required weight of 
coolant (and hence of missile) is obtained with a shape having a high 
pressure drag coefficient, that is to say, a blunt shape. On the other 
hand, 1t was found that 1f the missile were relatively heavy the least 
required weight of coolant, and hence of missile, is obtained with a 
shape having a low skin-friction drag coefficient, that is to say, a 
long slender shape. pa 


The second aspect of the heating problem treated was concerned 
with the rate of heat input, particularly with regard to thermal shell 


ORDA ta | AR 
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stresses resulting therefrom. It was seen that the maximum average 
heat-input rate and, hence, maximum average thermal stress could be 
decreased by using either a blunt or a slender missile, while missiles 
of intermediate slenderness were definitely to be avoided in this cor- 
nection. The region of highest local heat-transfer rate and, hence, 
probably greatest thermal stress was reasoned to be located at the 
forward tip of the missile in most cases. This was assumed to be 

the case and it was found that the magnitude of this stress was reduced 
by employing a shape having the largest permissible tip radius and 
over-all drag coefficient; that is to say, the blunt, high drag shape 
always appears to have the advantage in this respect. 


These results provide us with rather crude, but useful, bases for 
determining shapes of missiles entering the atmosphere which have mini- 
mized heat-transfer problems. If the over-all design considerations 
of payload, booster, et al, dictate that the re-entry missile be rela- 
tively heavy in the sense of this report, then it may be most desirable 
to make this missile long and slender, especially if the entry speed is 
very high (say 20,000 ft/sec or greater). Perhaps the slender conical 
shape is appropriate for such a missile. It seems clear, too, that the 
tip of this missile should be given the largest practicable nose radius 
in order to minimize the maximum local heat-transfer rate and hence 
maximum local shell stress problem. Even then it may be necessary to 
employ additional means to minimize the heat-transfer rate and, hence, 
thermal stress: encountered in this region (e.g., by sweat cooling). 


Let us now consider the case where the over-all design conditions 
dictate that the re-entry missile be relatively light in the sense of 
this report. This case is believed to be of more immediate importance 
than the one just considered since the lower sizes, weights, and 
entrance speeds to which it applies are more nearly in Line with those 
presently of interest. The relatively light re-entry missile will 
therefore be treated at greater length. 


A shape which should warrant attention for such missile applica- 
tion is the sphere, for it has the following advantages: 


1. It is a high drag shape and the frictional drag is only a few 
percent of the total drag. 


2. It has the maximum volume for a given surface area. 
3. The continuously curved surface is inherently stiff and strong. 


l. The large stagnation-point radius significantly assists in 
reducing the maximum thermal stress in the shell. 


Com 
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5. Aerodynamic forces are not sensitive to attitude and, hence, 
& sphere may need no stabilizing surfaces. 


6. Because of this insensitivity to attitude, a sphere may pur- 
posely be rotated slowly, and perhaps even randomly1“ 
during flight, in order to subject all surface elements to 
about the same amount of heating and thereby approach uni- 
Porm shell heating. 


On the other hand, the sphere, in common with other very high drag 
shapes may be unacceptable if: 


1. The low terminal speed permits effective countermeasures. 


2. The lower average speed of descent increases the wind drift 
error at the target. 


3. The magnitude of the maximum deceleration is greater than can 
be allowed. 


The first two of these disadvantages of the sphere might be minimized 
by protruding a flow-separation-inducing spike from the front of the 
sphere to reduce the drag coefficient to roughly half (see ref. 11). 
Stabilization would now be required but only to the extent required to 
counterbalance the moment produced by the spike. Special provision 
would have to be made for cooling the spike. 


These possible disadvantages of very high drag shapes may also be 
alleviated by another means, namely, using variable geometry arrange- 
ments. For example, an arrangement which suggests itself is a round- 
nosed shape with conical afterbody of low apex angle employing an 
extensible skirt at the base, as shown in figure 12. With the skirt 
Flared, the advantages of high drag are obtained during the entry 
phase of flight. As the air density increases with decreasing altitude, 
the skirt flare is decreased to vary the drag so as to produce the 
desired deceleration and speed history. If the deceleration is speci- 
fied in the equation of motion (see motion analysis), the required 
variation of drag coefficient with altitude can be calculated and, in 
turn, the heating characteristics can be obtained. 


14Note that if rotation is permitted, slow, random motion may be 
required in order to prevent Magnus forces from causing deviation 
of the flight path from the target. It should also be noted that 
at subsonic and low supersonic speeds gun-fired spheres, presumably 
not rotating, have shown rather large lateral motions in flight (see 
ref. 10). It is not known whether such behavior occurs at high 
supersonic speeds. 


dk 
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The examples considered, of course, are included only to demon- 
strate some of the means the designer has at hand to control and 
diminish the aerodynamic heating problem. For simplicity, this problem 
has been treated, for the most part, in a relative rather than absolute 
fashion. In any final design, there is, of course, no substitute for 
step-by-step or other more accurate calculation of both the motion and 
aerodynamic heating of a missile. 


Even from a qualitative point of view, a further word of caution 
must be given concerning the analysis of this paper. In particular, 
throughout, we have neglected effects of gaseous imperfections (such 
as dissociation) and shock-wave boundary-layer interaction on convec- 
tive heat transfer to a missile, and of radiative heat transfer to or 
from the missile. One would not anticipate that these phenomena would 
significantly alter the conclusions reached on the relative merits of 
slender and blunt shapes from the standpoint of heat transfer at 
entrance speeds at least up to about 10,000 feet per second. It can- 
not tacitly be assumed, however, that this will be the case at higher 
entrance speeds (see Appendix B). Accurate conclusions regarding the 
dependence of heat transfer on shape for missiles entering the atmos- 
phere at extremely high supersonic speeds must await the availability 
of more reliable data on the static and dynamic properties of air at 
the high temperatures and pressures that will be encountered. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. 28, 1953 
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APPENDIX A 
SYMBOLS 
A reference area for drag evaluation, Pt? 
B body factor, dimensionless 
(See eq. (51).) 
Cp drag coefficient, dimensionless 
Cr skin-friction coefficient based on conditions just outside 
the boundary layer, dimensionless 
Ce! equivalent skin-friction coefficient, dimensionless 
(See eq. (28).) 
TTL Lp 
Cp specific heat at constant pressure, ———>— 
slug “R 
Cy specific heat at constant volume, LER. 
slug OR 
F',F'',F'tt functions of PAy, dimensionless 
(See eqs. (18), (19), and (49).) 
g acceleration due to force of gravity (taxen as 32.2 rt.) 
sec 
ft-Lb 
b convective heat-transfer coefficient, 
ft? sec OR 
: T ft-lb 
H i | heat transferred per unit area, = 
Ft 
ft-lb 
k thermal conductivity, r 
sec ft? ([R/ft) 
m mass, slugs 
M Mach number, dimensionless 
Nu Nusselt number, dimensionless 
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Pr Prandtl number, dimensionless 
Q total heat transferred, ft-lb 
Re Reynolds number, dimensionless 
S surface area, ft” 
T temperature (ambient temperature of air at altitude y 
unless otherwise specified), OR 
t time, sec 
V velocity, It 
sec 
XI horizontal and vertical distance from impact point, ft 
2 
Z variable of integration, Z 
sec 
B constant in density - altitude relation, PL" 
(See ea. (2).) 
Y ratio of specific heat at constant pressure to specific 


heat at constant volume, Cp/Cy, dimenslonless 


A increment 
N distance within the shell measured normal to shell surface, fi 
9 angle of flight path with respect to horizontal, deg 
slugs 
u coefficient of absolute viscosity, —— S5 
ft sec 
slug 
p air density, 
TU” 
aj radius, ft 
y relative heat-transfer factor, dimensionless 


(See eqs. (56) and (57).) 
Subscripts 


o conditions at sea level (y = 0) 


i conditions at altitude ya(eg. (15)) 
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conditions at altitude y„(eq. (38)) 


conditions at altitude y,(eq- (#7)) 
conditions at entrance to earth's atmosphere 
local conditions 

recovery conditions 

stagnation conditions 

wall conditions 


conditions within the shell of the missile 
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APPENDIX B 


SIMPLIFYING ASSUMPTIONS IN THE CALCULATION OF 


AERODYNAMIC HEATING 


As noted in the main body of the report, the heating analysis is 
simplified by making the following assumptions: 


1. Convective heat transfer is of foremost importance; that is, 
radiative effects may be neglected. 


2. Effects of gaseous imperfections, in particular dissociation, 
may be neglected. 


3. Effects of shock-wave boundary-layer interaction may be 
neglected. 


l}. Reynolds! analogy is applicable. 
5. Prandtl number is unity. 


The restrictions imposed by these assumptions will now be considered 
in some detail. 


In assumption 1, two simplifications are involved; namely, (1) radia- 
tion from the surface of the body is neglected, and (2) radiation to the 
body from the high-temperature disturbed air between the shock wave and 
the surface is neglected. The first simplification may be justified on 
the premise that the maximum allowable surface temperature will be about 
the same for one body as compared with another, irrespective of shape, 
and, consequently, radiation away from the surface will be approximately 
the same. Hence, neglecting this form of heat transfer should not 
appreciably change the relative heating which is of principal interest 
in this paper. 


The second simplification of ignoring radiative heat transfer from 
the disturbed air to the body is not so easily treated. At ordinary 
Plight speeds this form of heat transfer is negligible since it is well 
established that at temperatures not too different from ambient tem- 
perature, air is both a poor radiator and a poor absorber. At the 
flight speeds of interest, temperatures in the tens of thousands of 
degrees Fahrenheit may be easily obtained in the disturbed air flow, 
especially about the heavier blunt bodies. At these temperatures it 
does not follow, a priori, that air is a poor radiator. Data on the 
properties of air at these temperatures are indeed meager. Hence, it 
is clear that calculations of radiative heat transfer from air under 
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these conditions must, at best, be qualitative. Nevertheless, several 
such calculations have been made, assuming for lack of better informa- 
tion that air behaves as a grey body radiator and that Wein's law may 
be used to relate the wave length at which the maximum amount of radia- 
tion is emitted to the temperature of the air (this assumption, in 
effect, enables low-temperature data on the emissivity of air to be 
used in calculating radiation at high temperatures). In these calcu- 
lations effects of dissociation in reducing the temperature of the 
disturbed air have also been neglected and hence from this standpoint, 
at least, conservative (i.e., too high) estimates of radiative heat 
trensfer should evolve. The results of these calculations indicate 

the Pollowing: (1) Radiative heat transfer from the disturbed air to 
the body is of negligible importance compared to convective heat trans- 
fer at entrance speeds in the neighborhood of, or less than, 10,000 feet 
per second; (2) Radiative heat transfer, in the case of relatively mas- 
sive blunt bodies, may have to be considered in heat-transfer calcula- 
tions at entrance speeds in the neighborhood of 20,000 feet per second; 
(3) Radiative heat transfer, in the case of relatively massive blunt - 
bodies, may be of comparable importance to convective heat transfer at 
entrance speeds in the neighborhood of 30,000 feet per second. From 
these results, we conclude, then, that the neglect of radiative heat 
transfer from the disturbed air to the body is probably permissible for 
all except, perhaps, very blunt and heavy shapes at entrance speeds up 
to 20,000 feet per second. However, this simplification may not be 
permissible, especially in the case of heavy blunt bodies entering the 
atmosphere at speeds in the neighborhood of, or greater than, 30,000 
feet per second. 


In assumption 2, the neglect of effects of gaseous imperfections, 
particularly dissociation, on convective heat transfer would appear to 
be permissible at entrance speeds up to and in the neighborhood of 
10,000 feet per second, since at such speeds the temperatures of the 
disturbed air are not high enough for these imperfections to become 
Significantly manifest. On the other hand, as the entrance speeds 
approach 20,000 feet per second, temperatures of the disturbed air may 
easily exceed 10,000° Rankine, in which case appreciable dissociation 
may be anticipated, inside the boundary layer for all bodies, and 
inside and outside the boundary layer in the case of blunt bodies. ‘The 
magnitude of these effects is at present in some doubt (see, e.g., the 
results of refs. 12 and 13.) Hence, for the present, the neglect of 
effects of gaseous imperfections on convective heat transfer is not 
demonstrably permissible at entrance speeds in the neighborhood of 
20,000 feet per second or greater. 


Tn assumption 3, it has been shown by Lees and Probstein (ref. 14), 
and more recently by Li and Nagamatsu (ref. 15), that shock-wave 
boundary-layer interaction may significantly increase laminar skin- 
friction coefficients on a flat plate at zero incidence and Mach 
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numbers in excess of about 10. Lees and Probstein found somewhat the 
opposite effect on heat-transfer rate in the case of weak interaction. 
It is not now known how this phenomenon depends upon body shape or type 
of boundary layer. However, it is reasonable to anticipate that there 
will be some effect, and certainly if the skin-friction coefficient is 
increased in order of magnitude at Mach numbers approaching 20, as 
indicated by the results of Li and Nagamatsu for strong interaction, 
then the phenomenon cannot be presumed negligible. Hence, we conclude 
that from this standpoint, also, the convective heat-transfer calcula- 
tions of this report may be in error at entrance speeds of the order 
of 20,000 feet per second or greater. 


The assumption that Reynolds' analogy may be used to relate skin- 
friction and heat-transfer coefficient does not, especially in the 
light of recent work by Rubesin (ref. 16), seem out of line with the 
purposes of this paper, at least at entrance speeds up to and in the 
neighborhgod of 10,000 feet per second. However, it does not follow, 

& priori, that this assumption remains valid at substantially higher 
entrance speeds, especially in view of the imperfect gas and shock-wave 
boundary-Layer-interaction effects already discussed. 


The assumption of Prandtl number equal to unity would also appear 
permissible for the analysis of relative heating of missiles at the 
lower entrance speeds considered here. However, in view of the ques- 
tionable effect (see again refs. 12 and 13) of dissociation on Prandtl 
number, it is not clear that this assumption is strictly valid at the 
intermediate and higher entrance speeds treated in this report. 


From these considerations it is concluded that the simplifying 
assumptions made in the main heat-transfer analysis of this paper will 
not significantly influence the results at entrance speeds in the 
neighborhood of or less than 10,000 feet per second. However, at 
entrance speeds in the neighborhood of and greater than 20,000 feet 
per second, these results must be viewed with skepticism. More accurate 
calculations of heat transfer at these speeds must, among other things, 
await more accurate determinations of both the static and dynamic prop- 
erties of air under these circumstances. 
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